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\^ . Arkani-Hamed, Dimopoulous, and Dvali have proposed a model of low-scale quantum gravity 

' featuring large extra dimensions. In this model, the exchange of Kaluza-Klein towers of gravitons 

I can enhance the production rate of electron and photon pairs at high invariant mass in proton- 

, antiproton collisions. The amount of enhancement is characterized by the parameter Ms, the 

, ^ ' fundamental Planck scale in the bulk extra dimensions. We have searched for this effect using 

' 100 pb~^ of diphoton data and 110 pb~^ of dielectron data collected with the Collider Detector 

' ' at Fermilab at ^/s = 1.8 TeV during the 1992-1996 run. In the absence of a signal, we place 95% 

confidence- level limits on Ms of 905 GeV/c^ and 826 GeV/c^, for the case of constructive and 
' destructive graviton interference respectively. 

PACS numbers: 04.50.+h, 04.80.Cc, 13.85.Qk, 13.85.Rm 

ly-^ \ Theories of low-scale quantum gravity featuring large extra spatial dimensions (LED) have attracted consid- 
^r^i • erable interest because of their possible observable consequences at existing and future colliders. In one such 
' scenario, proposed by Arkani-Hamed, Dimopoulos, and Dvali T'l, the fermions and gauge bosons of the Standard 
CO , Model (SM) are confined to the three ordinary spatial dimensions, which form the boundary ("the brane") of a 
' space with n compact spatial dimensions ("the bulk") in which gravitons alone can propagate. In this model, 
the Planck scale is lowered to the electroweak scale of 0(1 TeV), which is postulated to be the only fundamental 
scale in nature. The fundamental Planck scale in the extra dimensions {Ms), the characteristic size of the n 
^ ' extra dimensions (R) and the Planck scale on the brane (Mpi) are related via 
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Mil « M-+^R\ (1) 



a purely classical relationship calculated by applying the A + n dimensional Gauss's law. In this scenario, 
then, the weakness of gravity compared to the other SM interactions is explained by the suppression of the 
, gravitational field fiux by a factor proportional to the volume of the extra dimensions. 
; I ■ One important consequence for physics in the brane is that the discrete momentum modes of excitation 
I of the graviton transverse to the brane propagate in our three ordinary dimensions as different mass states. 
Analogously to the momentum states, the spacing between these mass states is proportional to 1/R. This 
collection of mass states forms what is known as a Kaluza-Klein (KK) tower of gravitons. The tower can in 
principle extend up to infinity, but there is a cutoff imposed by the fact the effective theory breaks down at 
scales above Ms- 

The existence of KK gravitons can be tested at colliders by searching for two different processes: real graviton 
emission and virtual graviton exchange. At leading order, virtual graviton exchange includes processes in which 
a virtual graviton is produced by the annihilation of two SM particles in the initial state, the graviton then 
propagates in the extra dimension and finally decays into SM particles that appear in the brane. Real graviton 
production occurs when a graviton is produced together with something else by the interaction of SM particles 
and escapes into the extra dimensions, leaving behind missing energy. 

In this paper, we search for LED by looking for the decay of KK gravitons into diphotons ,2i| or dielectrons 0. 
An analysis of these final states has previously been carried out by the D0 Collaboration A]. These channels 
are desirable because SM production of these pairs involves an electromagnetic coupling to the final state; thus 
production is suppressed with respect to dijet final states. However, gravitons couple democratically to stress- 
energy. Thus, LED production of dielectron and diphoton final states can more easily be distinguished from SM 
sources. The dominant Feynman diagrams that contribute to these processes in the SM are shown in Figure ^ 
When KK gravitons are included, new diagrams appear as shown in Figure |21 Since the LED contribution to 
SM pair production proceeds through a KK tower of graviton states with a closely spaced mass spectrum, the 
extra-dimensional signal does not appear as a single resonance, but rather as an enhancement of the production 
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FIG. 1: Feynman diagrams for SM diphoton and dielectron production at leading order. 




FIG. 2: Feynman diagrams for LED diphoton and dielectron production at leading order. 



cross section at high invariant mass where the SM contribution is rapidly falling and a large number of gravitons 
can be produced or, equivalently, more modes of the momentum in the bulk can be excited. 

The cross section for graviton production at leading order, as a function of the invariant mass of the electron 
or photon pair, takes the general form |^: 



^ = /sm(M) + VLEDflNTiM) + vlEofKKiM), (2) 

where M is the invariant mass of the photon or electron pair and fsM, fiNT, /kk are terms in the cross 
section with the dimensionful parameter tjled factored out. The term fsM on the right-hand side of Equation|21 
is the contribution from SM processes only. The second term, fiNT, proportional to the first power of tj^ed, is 
the result of the interference between the SM and the LED graphs. Finally, the last term, /kk, proportional to 
Vled^ is the contribution coming from the direct KK tower exchange. Graviton-mediated processes introduce 
the r/iBD-dependence in the equation, where tiled = ^/Mg is the parameter that characterizes LED: Ms 
represents the Planck scale in the bulk and should be of order 1 TeV and A is a dimensionless parameter of 
0(1) and its value is model-dependent. The different signs of A allow for different signs of the interference 
between SM and LED graphs. The shapes of the interference and direct KK cross sections are independent of 
Vled, which affects only the relative and absolute normalization, while the shapes themselves depend only on 
the parton distribution functions (PDF's) and the kinematics of the process. 

As an illustration, Figure 13 shows the diphoton invariant mass spectrum for the case Ms = 750 GeV/c^ with 
constructive interference between the SM and LED diagrams. The LED signal clearly stands out above the 
background at higher values of the invariant mass. The signal in the dielectron case is similar. We have searched 
for this signal in approximately 110 pb~^ of data collected in proton-antiproton collisions at ^/s = 1.8 TeV by 
the Collider Detector at Fermilab (CDF) experiment during the 1992-96 run (Run I). 

The Run I CDF detector has been described in detail in Reference 0. Here we present an overview of the 
detector components most relevant to this analysis. The central portion of the detector is cylindrically symmetric 
around the beam pipe. Forward/backward extensions together with the central section of the detector provide 
close to An coverage of the solid angle. Three tracking systems (SVX, VTX, CTC), situated inside a 1.4 T 
superconducting solenoid, allow measurements of charged particle momenta for pseudorapidity \ri\ < 1.1. In 
addition, the SVX allows measurement of secondary vertices and the VTX provides information on the location 
of the primary vertex along the beam axis. Located right outside the solenoid coil, the Central Preradiator 
(CPR), a multi-wire proportional drift chamber, measures tracks of charged particles that result from photon 
conversions in the inner detector. The electromagnetic (EM) and hadronic calorimeter systems are also located 
outside the solenoid coil, past the CPR. They consist of sampling devices with alternating layers of sensitive 
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FIG. 3: An example of LED {Ms ~ 750 GeV/c^) contributions is shown along with the SM prediction in the diphoton 
invariant mass distribution. The two curves for LED are for constructive (A — —1) and destructive (A — +1) interference 
with the SM processes. 

and absorbing materials, and are segmented in r/, towers with projective geometry. The calorimeter is divided 
into three components covering different regions of pseudorapidity: central (| ry |<1.1), end-plug (1.1<| 1] |<2.4) 
and forward (2.4<| 77 |<4.2). Only the central and end-plug regions are used in this analysis. The central 
strip chambers (CES) are embedded into the central EM calorimeter section, and provide a high-resolution 
measurement of the lateral shower profile in z and r — at the point of its maximum development. Finally, 
muon chambers are located outside the calorimeter to detect particles that escape the inner layers of the detector. 

The diphoton sample consists of two distinct subsets, corresponding to two different event topologies. One 
of the subsets, the Central-Central (CC) diphoton sample, consists of events with both photons in the central 
region of the detector, while the other, the Central-Plug (CP) diphoton sample, has events with at least one 
central electromagnetic cluster and at least one cluster in the end-plug electromagnetic calorimeter. The CC 
diphoton sample has been described in Reference 0. Photons in this sample are required to have transverse 
energy (Et) > 22 GeV, to pass a calorimeter isolation cut, to have no nearby charged tracks, and to have a 
CES shower profile consistent with that of a photon. The identification efficiency of the CC diphoton selection 
is 0.626±0.017, and the integrated luminosity of this sample is 100 pb"-'^. A total of 287 CC events pass these 
cuts. 

Two trigger paths are used to select the CP diphoton sample. The first path requires at least one isolated 
central EM cluster with Et > 23 GeV, while the second path requires an Et > 50 GeV central EM cluster 
with no calorimeter isolation requirement. The central cluster must pass the the same offline selection applied 
to the CC diphoton sample but with an Et > 25 GeV to ensure a high trigger efficiency. The plug EM cluster 
is required to have Et > 22 GeV. CTC tracking does not extend to the plug region of the detector; thus to 
reject charged particles we require that the fraction of hit wires in the VTX in a road centered on the cluster 
direction is smaller than 0.4. We require the Et in the annulus between 0.4 and 0.7 in ij-cj) space around the 
cluster to be no larger than 4 GeV. This requirement allows us to select isolated clusters while avoiding possible 
inefficiency due to additional leakage into the cone immediately surrounding the cluster that occurs at larger 
rapidities. The efficiency of this cut is calculated with electrons from a Z — > e+e~ control sample, and is found 
to be 89%, 10% more than the canonical isolation cut based on a cone of 0.4. The identification efficiency of the 
CP diphoton selection is 0.550±0.065. In our integrated luminosity of 87 pb^^, 192 CP events pass our cuts. 

Dielectron events used in this analysis are required to have at least one "tight" electron candidate in the 
central calorimeter, while the second "loose" candidate can be in either the central or plug calorimeter. The 
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details of the selection requirements are described in References [3|, |8j . All events must pass a central electron 
trigger that requires an EM cluster with > 22 GeV together with a track with px > 13 GeV/c. Offline, the 
tight electron must pass an Et cut of 25 GeV and have the ratio of energy to momentum {E/p)< 4. For tracks 
with pt > 50 GeV/c, no E/p cut is imposed. We require good agreement between the EM shower location and 
the extrapolated track position, and require the ratio of hadronic to electromagnetic energy in the cluster to be 
low. Finally, wc require the electron to be isolated. The loose EM cluster is also required to have Et > 25 GeV, 
but central candidates may pass more relaxed isolation and hadronic energy cuts. Plug electron candidates are 
required to have a low hadronic energy fraction, to be isolated, and to have a lateral shower profile consistent 
with that of testbeam electrons. The efficiency of the electron identification cuts is measured using Z data and 
is found to be 0.916 ±0.010 for plug electrons and 0.961 ±0.024 for central electrons. In an integrated luminosity 
of 110 ph-\ we identify 3319 CC and 3825 CP candidates. 

Both diphoton and dielectron events are generated using a leading-order matrix-element Monte Carlo (MC) 
program jOj]. The simulation includes the matrix elements for direct SM diphoton production with qq or gg 
in the initial state, production via KK exchange from a qq or gg initial state, as well as interference between 
the SM and LED diphoton production processes with qq in the initial state. However, it does not include 
the interference between the g^-initiated SM and LED initial states, which is expected to be small jlOjj. A 
conservative estimate of this contribution is incorporated into the systematic uncertainties. The simulation 
also generates SM qq — > e^e~ and qq or gg — > e^e~ events mediated by a KK tower of gravitons. The 
PYTHIA MC program is used to model fragmentation, parton showers, and the underlying event with 
CTEQ5M PDF's The generated events are passed through a fast simulation of the CDF detector. We 

select simulated events using the same cuts applied to the data. This procedure yields a prediction of 96 ± 31 
prompt SM diphoton events in the CC sample and 76 ± 14 events in the CP sample. In the dielectron sample, 
the prediction is 3463 ± 223 CC and 3883 ± 292 CP prompt SM dielectron events. 

However, the bulk of the diphoton sample consists of events in which one or both electromagnetic clusters 
originate from jets. Hard jet fragmentation to leading neutral mesons (7r*',7y) that subsequently decay into 
multiple photons can fake the prompt photon signature. In the CC case, the fake background is evaluated by 
employing statistical methods 12] that use CES and CPR pulse shape and height information to discriminate 
between prompt photons and photons originating from neutral meson decays. These techniques yield a fake es- 
timate of 64±11±19%, which corresponds to 183±56±32 events. The mass spectrum of the CC fake background 
is determined from the invariant mass distribution of events in a control sample of non- isolated diphotons . 

In the CP case, we rely on the above techniques to determine the probability that the central cluster is a jet. 
However, these methods cannot be used for plug clusters. Instead, we employ dijet and non-isolated diphoton 
control samples to estimate the fake rate of events in which the central cluster is a prompt photon while the 
plug cluster is a jet 0. The total number of fake events in the CP sample is 132±31. The mass spectrum of 
the CP fake background is determined by simulating SM 7+jet processes using the PYTHIA MC program and 
the CDF detector simulation. These events provide the mass spectrum for 7±fake events. The mass spectrum 
for events where both clusters are fakes is obtained by multiplying the MC 7±jet shape by the ratio between 
the 7r°7r" and 77r" cross sections calculated at leading order 13] . Because tt^'s are the main contribution to the 
fakes, we take this rescaled shape as the shape of the fake±fake contribution. 

In summary, in the diphoton channel the prediction of fakes plus SM 77 is 280±66 events in the CC sample 
and 208±34 events in the CP sample, in good agreement with the observed 287 events in the CC case and 192 
events in the CP case. 

In the dielectron sample, the electrons in the final state are produced by other SM processes in addition 
to Drell-Yan. The dominant contribution is from QCD dijet events in which each jet contains a real or fake 
electron that survives the electron identification cuts. This mass distribution and size of this background is 
evaluated using samples of non-isolated dielectrons. We obtain 10±6 background events in the CC case and 
224±17 background events in the CP case. 

In summary, in the dielectron channel the number of events predicted, adding together background and SM 
contributions, is 3473 ± 223 CC and 4107 ± 293 CP dielectron events, in agreement with the observed 3319 
events in the CC case and 3825 in the CP case. 

We find that the data are in good agreement with the SM plus fakes prediction both with respect to the 
rate and invariant mass distribution. Figure 01 shows the data overlaid with the background predictions for 
diphoton and dielectron events. We therefore proceed to set a limit on LED by determining the lower bound 
on the parameter Ms that characterizes the strength of the graviton coupling. We do this by performing a 
likelihood fit to the diphoton and dielectron invariant mass distributions. We use two very similar likelihood 
functions, one for the diphoton and one for the dielectron case, which are multiplied together to obtain the 
overall likelihood function C. The form of £ is similar to that employed in CDF's measurement of the top quark 
mass 0] . In the likelihood function, a Gaussian constraint is applied to the number of SM and fake dielectron 
and diphoton events in the CC and CP samples. In addition, a Poisson constraint is applied to the total number 
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of dielectron and diphoton events. Finally, a shape term weights the invariant mass of each event according 
to the consistency of its invariant mass with the SM, backgromid, and LED shapes. We minimize — ln£ to 
obtain the best-fit value of tjled = A/M| for the cases A = ±1, corresponding to constructive and destructive 
interference of the SM and KK terms. These values are consistent with zero, so we extract 95% confidence level 
(C.L.) limits. For the diphoton data alone, the limits on vled for constructive and destructive interference 
respectively are "qLED^^-i^ TcV^* and riLED=2.18 TeV^"'. For dielectrons, the limits are tiled— "^-QQ TeV^''' 
(constructive) and tiled— '^■SG TeV^"* (destructive). 

Systematic uncertainties in this analysis arise from effects that can modify the shapes of the SM, fake, and 
LED invariant mass distributions. In both the electron and photon analyses, the dominant effects are the choice 
of PDF's, variations in initial state radiation (ISR), and the shape of the CP fake background. Furthermore, 
the interference between SM and LED gg initial states is not included in the diphoton MC. To quantify these 
systematics, we generate shifted MC invariant mass distributions with each systematic effect separately varied. 
For example, we vary the choice of PDF's (CTEQ5HJ, MRST(h-g), MRST(larger d/u ratio), MRST(low-a^), 
MRST((7rf), instead of the default CTEQ5M ^^'^ turn off the effects of ISR. To estimate the uncertainty 

from the neglected gg interference term, we note that the cross section for the gg — > 77 process is about 30% 
of the qq process, so the gg interference term should be roughly (0.3)^ w 0.1 of the qq one. We vary the 
qq interference portion of the templates by twice this amount, ±20%. The systematic uncertainty on the CP 
background is estimated by using the second most populated region of the non-isolated dielectron control sample 
to determine the background shape, and computing the change in fitted ti^ed versus the default (most populated 
region). For each systematic effect, many pseudo-experiments are generated at different values of AIs from the 
shifted distributions, and the resulting invariant mass distributions are fitted to our default templates. For each 
systematic effect, we determine the average shift between the input value of tiled and the value extracted from 
the fit. The resulting shifts for all systematic effects are added in quadrature. The overall shift is used as the 
width of a gaussian that is convolved with the raw likelihood function from the fit to the data. Remaining 
uncertainties on efficiencies, acceptance, luminosity, and background normalization are already included in the 
width of the gaussian constraints in the likelihood function. 

In the diphoton analysis, the 95% C.L. limits with systematics included are 77l£;c = 1.53 TeV^* (constructive 
interference) and 77^^0=2.48 TeV^"' (destructive). In the dielectron case, we find riLED='^-7Q TeV^"' (construc- 
tive) and TiLED=2-87 TeV~^ (destructive). 

Finally, we obtain the strongest limits by combining the diphoton and dielectron samples. A new fit is per- 
formed by multiplying the likelihood functions for each sample and minimizing it with respect to the common 
parameter tiled- The systematic uncertainties are included into the combined limit in the following way. The 
full systematic uncertainties for the two search channels, Atiled eeiVLED) and Atiled 'ryiVLED), are broken 
into their component pieces, Sv£EE,ee, Svi!h%ee, ^dX- Sv£EE,y^, ^^1%%^^ and StiH^S^W^ Here SvH'^^;^ 
represents the sum in quadrature of the dielectron systematics with the diphoton systematics, i.e. background 
and LED acceptance. Similarly, STiYEo^-t represents the sum in quadrature of the diphoton systematics un- 
correlated with the dielectron systematics, i.e. SM-t-LED interference cross section and CP background. Each 
likelihood function, diphoton and dielectron, is smeared by taking the corresponding, unsmeared likelihood 
function, and smearing each point by throwing random numbers Ui, i=l,...,4, and smearing a sample of values 
of tiled iu the following way: 

VLED.ee = VLED.ee + (^I^VlED ,ee{VLED ,ee) + (^l^vl^ED .eeiVLED ,ee) + Cr3STll"£l^''^^{TlLED,ee) (3) 
VLEDni = VLEDrn + '^i^'nLED,-i'y{VLED,-i'y) + Cf25VL^ED,'i-l{VLED,-i-t) + (4) 

This ensures that the uncorrelated uncertainties are smeared independently while the correlated PDF's and 
ISR uncertainties are smeared together. 

After the smearing procedure is performed, the diphoton and dielectron smeared likelihood functions are 
multiplied together and the same method described above is used to calculate the combined dielectron-diphoton 
95% C.L. limit on vled- We find tiled— ^-^^ TeV~^ (constructive) and tiled='^-^5 TeV~^ (destructive). 

The limits on vled are converted into limits on Ms by using the equation vled = A/M|. In Hewett 
convention A = —1 for constructive interference and A = -1-1 for destructive interference. Other popu- 
lar conventions are those of Giudice, Rattazzi, Wells (GRW) ^17], and Han, Lykken, Zhang (HLZ) [l^- To 
translate results from Hewett to GRW convention one simply multiplies M5 (Hewett) by \/tt/2 (constructive 
interference only). In HLZ convention the dependence on the number of extra dimensions n is calculated and 
it is incorporated into A. For n > 2: 




Ms (HLZ) = ( - A^s (Hewett). (5) 
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FIG. 4: Invariant mass distributions for the diphoton (top) and dielectron (bottom) data overlaid with the expected 
invariant mass distributions from SM and fake sources, as well as the LED contribution for M5(Hewett) = 900 GeV/c^ . 



The dielectron and diphoton combined limits in these conventions are summarized in Table ^ The limits 
quoted so far assume a ii'-factor of 1.0 for the LED process. Limits can be improved somewhat by assuming 
if = 1.3 to account for next-to-leading order effects for the SM diphoton and dielectron production as well as 
for the corresponding graviton mediated processes, as was done in Reference jj]. These limits are also shown in 
Tabled 
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TABLE I: Diphoton and dielectron limits on Ms in Hewett, GRW, and HLZ conventions. The final two rows show the 
combined limits with and without a i^-factor. 
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